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^ . I. INTRODUCTION 

o , 

In the analysis of General Relativity as a Hamiltonian system [1] one chooses a time function and considers the 
foliation of the spacetime by the slices of constant time. Two natural geometrical quantities arise on such spacelike 
three slices. One is the intrinsic three metric, usually g ab , and the other is the extrinsic curvature K ab , the derivative 
of g ab along the normal to the slice. They are related by the constraints, which in a vacuum spacetime read 

ft(3) _ K ab K ab + (trK) 2 = (1) 

VaR ab _ g ab VailR = (2) 



^1" 
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where is the three scalar curvature. Given the initial data one chooses, essentially arbitrarily, the lapse, N, and 
the shift, N l , which determine the magnitude and direction of the unit time vector relative to the normal to the slice. 

One now can write the evolution equations for the intrinsic metric and extrinsic curvature in vacuum, e.g. [2] (the 
reader should be warned that we follow Wald [3] in our definition of the extrinsic curvature, not [1], positive K means 
■ inceasing volume to the future) 

m : 

d t gab = 2NK ab + N a:b + N h . a (3) 
d t K ab = N. ab - N (R ab - 2K d a K bd + K ab trK) + K ab . c N c + K ac N? b + K cb m a . (4) 

a -1 

Let us stress that we are using the convention of signs that gives trK = +n a a where n a is the timelike unit normal to 
&)[)■ the slice and dt^fg = y/g (NtrK + N.fA. In other words, positive trK means expansion. It is often useful to specify 
the foliation, and thus the time, by placing a condition on the extrinsic curvature. A very popular choice is to demand 
• ■ *h , that the trace of the extrinsic curvature be constant on each slice ('CMC slicing'). 

In this paper we investigate the CMC slices of the extended Schwarzschild solution. The manifold consists of four 
segments, each of which can be covered by the standard Schwarzschild metric 

ds 2 = -(l-^j dt 2 + + r 2 (d9 2 + sin 2 9d<f> 2 ) (5) 

where t is the 'static' killing vector and r is the 'areal' or 'Schwarzschild' radius. In the left and right zones t is timelike 
and r > 2m is spacelike. In the bottom zone r is timelike and runs forward from the past singularity at r = to 
r = 2m. In the top zone r is also timelike and runs forward from r = 2m to the future singularity, also at r — 0. We 
do not seek the most general CMC slices. We are looking for those CMC slices which inherit the underlying spherical 
symmetry of the given spacetime. 

There are two complementary ways of analysing this problem. One way is to assume that one is given initial data 
(the intrinsic metric and extrinsic curvature) both parts of which have the desired symmetry. It turns out that one 
can explicitly solve the constraints. From the momentum constraint, Eq.(2), it is clear that the extrinsic curvature 
must be just a sum of the trace term and a part which is both trace and divergence free ('transverse-tracefree', 'TT'). 
There exists a unique spherically symmetric TT tensor. Therefore the extrinsic curvature can be written down with 
just two free parameters. On substituting into the Hamiltonian constraint, Eq.(l), one discovers that this also can be 
solved explicitly. 

The alternative approach is to take the given spacetime and make a coordinate transformation in the (t, r) plane 
only, given by t = h(r), leaving the rest untouched. h(r) is called the height function. One now imposes the condition 
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that the t' = slice be CMC. This gives a second order equation for the height function which can be integrated 
explicitly once. This is enough to evaluate the intrinsic metric and extrinsic curvature of the slice and, of course, they 
agree with the expressions obtained using the first approach. 

One then can work out how these slices fit into the given spacetime and construct interesting CMC foliations. One 
ends up with a first order equation for height function which cannot be integrated explicitly. Nevertheless one can 
make qualitative statements about the location of the slices. In this article we focus on a particular class of slicings 
where we fix the value of tiK and vary the parameter defining the amount of the TT component in the extrinsic 
curvature. For small values of the parameter we have two foliations, one which runs from one infinity to the other, 
and one which emerges and returns to one of the singularities. As the parameter increases, the leaves of the foliations 
approach one another, and at a critical value of the parameter, they touch. For values of the parameter greater than 
the critical one the nature of the CMC slices change. They all now run from infinity into the singularity. 

In this article we focus on the behaviour of the slices as they approach the critical value. We find the classic 'collapse 
of the lapse' phenomenon. Further, by looking carefully at the first order equation for the height function we obtain 
an explicit expression for how the lapse decays near criticality. This article draws very heavily on the analysis given 
in [6] of the collapse of the lapse for maximal slices of Schwarzschild. 

II. EXTRINSIC AND INTRINSIC GEOMETRY OF CMC SLICES 

We can generate an essentially general spherically symmetric slicing of the Schwarzschild solution by making a 
coordinate transformation t' = t'(t, r),r' — r'(t, r). This will give us a spacetime metric of the form 

ds 2 = -gt>t>dt' 2 + N r ,dt'dr' + g z > z ,dr' 2 + r 2 (d9 2 + sm 2 (6)d(f> 2 ) , (6) 

We have that gt>t',N r >, and g r ' r i are functions of t 1 and r' . The coefficient, r 2 , in front of the two-metric is the 
original, Schwarzschild, coordinate r 2 but can be viewed as a function of t' and r' as well. We can make a further 
coordinate transformation of the form r" = r"(t',r'), leaving t' unchanged. This has the effect of changing the r' 
coordinate within each slice but leaving the slicing unchanged. This kind of transformation can be used to arrange 
that Vi' • Vr" = 0. This is equivalent to dragging the r' coordinate along the normal to the slicing and thus sets the 
shift to zero. This will give us a spacetime metric 

ds 2 = -N 2 dt' 2 + g r „ r „dr" 2 + r 2 (d0 2 + sm 2 (0)d^ 2 ) . (7) 

On any one given slice we can arrange that r" = r, the original Schwarzschild coordinate. However, when one tries to 
propogate this condition one discovers that it is not compatible with vanishing shift. Therefore one can have a metric 
of the form (6) with r' = r or a metric of the form (7) with zero shift but not both. Another choice would be to set 
g r i r i = 1, i.e., to choose the r' coordinate as the proper distance along the slice. Again this is not compatible with 
vanishing shift. One advantage that the 'proper distance' coordinate choice has over the V = z' gauge is that, so long 
as the slice remains spacelike, the proper distance gauge always remains regular while the V = r' choice may well 
have coordinate singularities. However, in this article we will largely stick to the metric form (7) and ignore questions 
such as the 'best' choice of spatial coordinate. 

To simplify the notation we will write the line element as 

ds 2 = ~N 2 dt 2 + adr 2 + R 2 (d0 2 + sm 2 (0)d<j) 2 ), (8) 

where the written (t, r) are NOT the original (t, r) while R IS the original r. The geometry of t =const slices is 
encoded in two places. One is the dependence of a on r and the other is the relationship between R and r. This 
second piece is contained in the mean curvature of the surfaces of constant r as embedded two-surfaces in the spatial 
three-geometry, 

2 dR 
^ yfaR dr 

The only nonzero (three) extrinsic curvature components with mixed-case indices are 

K ; = ^1, K $ = k* = M = I (tr # _ K ). (9) 
r 2aiV' 6 * RN 2 V r> K ' 

These can be viewed as evolution equations for a and R. The evolution equations for the extrinsic curvature can most 
compactly be written as 
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d t txK = V«V 4 iV - K)K\N (10) 



d t {tiK - K r r ) = NR^ r - NK{K) + NKltiK + (11) 

V d 



The form of the three-dimensional Ricci curvature component R^ r is 

%(pR) 



R?» = (12) 



while the three-dimensional scalar curvature R^ is 



(3)= 2d r (pR) ( P R) 2 2 
2i? 2 R 2 ' 



(13) 



It turns out that the momentum constraint can be written as 



d r {K r r - tvK) = - ^pK r r + iptrtf, (14) 



and the Hamiltonian constraint is 



^d r (pR) = -\(Kf \p 2 + ± + + i(tri,) 2 . (15) 

We are interested in finding surfaces which have trK = constant, where trK = (1/N)dt h\(y/aR 2 ) is the fractional 
rate of change of a coordinate volume during the temporal evolution. Assume that K = trK is constant on a fixed 
Cauchy hypersurface. Then the momentum constraint (14) is easily solved by 

Kr = J + W , (16) 

where C is again a constant on the chosen Cauchy slice. The other components of the extrinsic curvature are 

*»-*J-f-»- < 17 » 

This can be recognised as a combination of the trace term plus the unique spherically symmetric TT tensor, the terms 
with coefficient C. Therefore CMC slices of Schwarzschild are completely defined by the two parameters K and C. 
The only residual freedom is the ability to drag any surface along the Killing vector without disturbing either the 
intrinsic or extrinsic geometry. 

The insertion of (16, 17) into the Hamiltonian constraint leads, after some minor manipulation, to the equation 



d r 

Eq.(18) is solved by 



R (nR) 2 ff ° 2 K \ 3 



0. (18) 



<*>*-'H + (™-£))- 

Here (3 is essentially the integration constant, modified by completing the square of K and C related terms. It is easy 
to show that (3 — 2m. If 'r' is replaced by the 'areal radius, 'R' then one finds a = A/{pR) 2 . Notice also that the 
three-dimensional line elements reads in such a case 

ds 2 3) = j^dR 2 + R 2 (d8 2 + sin 2 (0)# 2 ). (20) 
ypri) 
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III. THE CYLINDRICAL CMC SLICES OF THE SCHWARZSCHILD SPACETIME 



In the upper and lower quadrants of the Schwarzschild spacetime the Killing vector is spacelike and runs along the 
r = constant surfaces. Since everything is constant along the Killing vector, the trace of the trace of the extrinsic 
curvature is preserved along these cylindrical surfaces. Therefore each r = constant surface is a CMC slice. 

The trace of the extrinsic curvature (in the upper quadrant) is given by 

K= 2r - 3m (21) 
V2mr 3 — r 4 

This is large and positive near r = 2m, zero at r = 3m/2 and becomes large and negative as r becomes small. This 
transforms into 

r 4 _ 2mr 3 + ^r-irnf = Q ^ 

This is a quartic equation with two real roots. One lies between r = 2m and r = 3m/2 and the other between 
r = 3m/2 and r — 0. This is clear by looking for the extrema of the quartic. To find these we just differentiate to get 

o r 3m 2 

4r 3 - 6mr 2 + 4—^ — = 2(2r - 3m) (r + = 0. (23) 

Therefore it has only one minima (at r — 3m/ 2) and the quartic is negative there. Hence it has two real roots, one 
on each side of r = 3m/2. On substituting back into Eq.(21) it is clear that the solution of Eq.(22) with r > 3m/2 
has K > and the solution with r < 3m/ 2 has K < 0. 

In the lower quadrant things are somewhat different. The trace of the extrinsic curvature is now given by 

K—S^-,. (24) 
v2mr J — 

This is because, in the upper quadrant, the future is in the direction of decreasing r while in the lower quadrant the 
future is in the direction of increasing r. This now is large and positive near r = 0, zero at r = 3m/ 2 and becomes 
large and negative as r approaches 2m. We get the same quartic, Eq.(22), with the same roots, but now with order 
reversed. The root which is less than 3m/ 2 corresponds to K > 0, while the other root has K < 0. 

Given r, we can work out, from Eq.(21), the value of the trace of the extrinsic curvature. We can, in fact, work out 
the entire extrinsic curvature and evaluate the constant C associated with these cylindrical CMC slices. In the upper 
quadrant we get 

3\f2mr z — r 4 ^ ^ 

Therefore C > on each of these slices. 

In the bottom quadrant, the extrinsic curvature picks up a minus sign. Therefore for the cylindrical CMC slices we 
get 

K= 2r-3m _ g = r 4 - 3mr 3 

\j2mr z — r 4 3\j2mr z — r 4 ' 



and so C is negative in the lower quadrant. 



IV. THE EMBEDDING OF CMC SLICES IN THE SCHWARZSCHILD SPACETIME 

In addition to the cylindrical CMC slices described in Section III, there are many other spherically symmetric ones. 
In this section, we will discuss how they run through the spacetime. 

The first comprehensive analysis of CMC slices in Schwarzschild appeared in [5]. The analysis given here closely 
follows the analysis of the related problem of maximal slices {tvK = 0) in the Schwarzschild solution given in [6] . Let 
us start off with the Schwarzschild metric 

ds 2 = - (l - ^) dt 2 + + r 2 (d9 2 + sin 2 Odcf) , (27) 
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and look at the slice given by t — h(r), where h(r), for obvious reasons, is called the height function. One way of 
understanding the geometry of this slice is to make the following coordinate transformation: 



t = t- h(r) 
f = r 
6=6 



t = t + h(f) 
r = f 



where the t = surface is the slice we are interested in. The transformed metric becomes 

/ -(l-M), - h >(i-^) o \ 

h' 2 (1 - 



9n» 



V 



-fc'(l-M), (l-M)- 













r z sin 



e) 



r v = 



/-(l-2p)- 1 +/i' 2 (l-22i) ) \ 

-/i'(i-^f), (i- 2 ^) 
o o^o 
V o oo -j-fen 



where ft' = dh/dr. The intrinsic metric is given by 

2m\ _1 



dr 2 + r 2 (dfl 2 + sin 2 Odcf) . 



(28) 



The lapse TV of this slicing is given by 



N 



2m\ 1 



r ) 



2m 
r 



the shift N a by 



N n = 



^).o.o 



and the future-pointing unit normal by 



(29) 



(30) 



(31) 



Given any three slice in the four manifold, we can drag it along by the killing vector. This will give a slicing where 
the time translation is along the killing vector. It is this slicing that is generated by the coordinate transformation 
above. Therefore the N and N a defined by Eqns.(29) and (30) arc nothing more than the projections of the killing 
vector perpendicular to and onto the given slice. Of course, the slicing given by dragging along the killing vector 
cannot form a foliation because the killing vector has a fixed point on the bifurcation sphere. 
The mean curvature of the t = slice is given by 



2 h >H_2m 
V r 



If K is a constant this can be integrated to give 



Kr C 



(I- 2 ?)" 1 -h>* (1-2?) 



(32) 



(33) 
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where C is a constant of integration. In turn, this can be manipulated to give 



i-^y - h *(i-*») = i „ (34) 

r J V r J (i_M) + (*r_ cf 

and hence 

Kr C 

ti = 3 g (35) 

If one could integrate this one more time and find h(r) in closed form, we would have a complete description of the 
slices. We cannot do so. Nevertheless, we can extract a significant amount of information from Eq.(35) as it stands. 

First, from the expressions Eq.(28), Eq.(29), and Eq.(30) it is clear that the intrinsic metric, the lapse, and the 
shift depend only on h n . Thus we get 



pr / / 2t7j\ / Kr C 



1 i T= \/( : -) ' • • - = - 



Finally, we can find the extrinsic curvature of the slice by using 



2NK ab = ^ - N a , h - 7V fc;a , (37) 

and we recover Eqns. (16) - (20). 

When we look at Eq.(35), it is clear that the right hand side does not decay for large r and thus the integral must 
diverge as we approach infinity. This is not surprising as we expect the CMC slices to go to null infinity. 

This agrees with the behaviour of the spherical CMC slices in Minkowski space. Consider the mass hyperboloid 
t 2 — r 2 = 9/K 2 in flat spacetime. If we choose the one which goes to future null infinity then the future-pointing 
timclike normal is n a — (t, r)j\Jt 2 — r 2 , where we have to choose the positive root of \Jt 2 — r 2 . We then find 
XxK = n a a = 3/\/t 2 — r 2 = \K\ > 0. To find out where on null infinity the slice ends up, we need to introduce 

null coordinates v = (t + r)/2,u = (t — r)/2. Using t = r^Jl + 9/r 2 K 2 w r + 9/2rK 2 it is clear that as v — > oo 
u w 9/rK 2 — > 0. If we time translate it to (t — to) 2 = r 2 + '9/K 2 we find u — > to- Therefore it slides up null infinity. 

If we look at Eq.(35) for large r we see that hi ~ rj \j9jK 2 + r 2 . The integral of this is h w ^/9jK 2 + r 2 which 
is in complete agreement with the flat space expression. Therefore the slices remain spacelike but go null infinity as 
r — > oo. Further, if K > the slices all go to future null infinity whereas if K < the slices go to past null infinity. 

One place we can find interesting information, without solving for h, is by looking at the expression for the mean 
curvature of the spherical two-surfaces. In particular, we know that 

p 2 r 2 2m\ t f Kr C x 2 
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= (1- — ) + (—--) >0- (38) 



Therefore the polynomial on the right hand side of Eq.(38) must be nonnegative. Further, we know that the zeros 
of the polynomial are the points where p = and therefore are the extrema of the area of the round two-spheres as 
embedded surfaces in the three-slice. 

Let us first fix some K > and see what happens as we vary C. (The cases where K < is remarkably similar.) 
First consider the case where C = 0. This is the so-called umbilical slice, where K a b oc g a b- In this case the 
polynomial reduces to 1 — 2m jr + K 2 r 2 /9. This is a cubic equation with only one real root, call it r u . Outside r = r u 
the polynomial is positive, inside it, it is negative. It is clear that r u < 2m and that r u = 2m iff K = 0. Therefore 
we know that the umbilical slice with K > starts out at future null infinity, comes in to a minimum at r = r u 
and then passes out to the other future null infinity. The obvious question is whether this occurs above or below the 
bifurcation sphere. 

To settle this we need to look at the optical scalars [7-9] 




w + = 2\i[ I- — )- (— - -) F2( — --). .39) 

and 



G 



(40) 



These are essentially the null expansions in the outgoing-future and outgoing-past directions respectively. They are 
both positive in Minkowski space and in the exterior regimes of the Schwarzschild solution. Since the product satisfies 
o! + w_/4 = f — 2m/r one or other becomes negative in the interior quadrants of the Schwarzschild solution. It turns 
out that the upper quadrant satisfies lj + < 0,a->_ > while in the lower quadrant we have that < 0,c<j + > 0. It is 
clear that at r = 2m, when K > and C = 0, ll>- goes negative while uj + remains positive. Therefore the umbilical 
slice (with K > 0) must pass through the lower quadrant. Therefore it starts at future null infinity, comes down so as 
to cross the t — axis, passes through the Schwarzschild throat below the bifurcation point to some minimum radius 
r u and then rises up again to the other future null infinity. 

Let us now hold K > fixed but change C so as to be slightly larger than zero. Now the polynomial becomes sixth 
order with two roots which we call r ms (ms = max. -small) and r m i (ml = min. -large). Near r — the dominant term 
is the positive term C 2 /r 4 so the polynomial starts off large and positive while then the next term is the negative 
— 2m/r which pulls it negative at r = r rns with r ms w y/ C' 2 /2m. We know that the polynomial must become positive 
before r — 2m and the K 2 r 2 /9 term does just that at r m i with r m i w r u . If C > then the effect of the C term is 
to diminish the effectiveness of the K 2 term so we get that r m i > r u while if C < we get r m ; < r u . Therefore for 
K > and C > (but small) we have two different regimes in which the polynomial is positive. One is for small 
r which represents a CMC slice which starts at r = 0, expands out to r ms , which is the maximum area, and then 
contracts again back to r = 0. 

When we look at the null expansions it is clear that when C > for small r we have that oj+ < 0, u>- > so 
it must be in the upper quadrant. Hence when C > the small r slice comes out of and goes back to the future 
singularity while the other slice runs from future null infinity to future null infinity and passes through the center at 
a slightly larger radius than the umbilical one. Thus it is to the future of the umbilical slice and crosses closer to 
the bifurcation sphere. As C increases away from zero we continue to have two CMC slices, one which comes from 
the future singularity out to some small radius r ms w y/ C' 2 /2m and the other which goes from future null infinity 
to future null infinity but will be slightly to the future of the umbilical slice. We find r TO ; monotically increases as 
C does until C = 8Km 3 /3. For this value of C, it is easy to show that r m i = 2m so that this CMC slice will pass 
through the bifurcation point. 

Increasing C acts like a time translation near infinity. From what happens in Minkowski space, we expect the slice 
to slide up along null infinity. 

As C increases even further r m i will start to decrease again while the CMC slice continues to move forward in time 
and passes through the throat above the bifurcation point. As C increases we find that r ms increases so that the CMC 
slicing that begins and ends at the future singularity moves backwards in time. The mimimum of the polynomial rises 
up and the two roots, r TO ; and r msi will approach each other as C approaches the critical value C = C+. For this 
value of C the polynomial is everywhere positive except at one point. This will be at a radius we call This will 
satisfy Rf > 3m/2. Rf is nothing more than the larger of the two roots of Eq.(22) and C+ is the value of C given 
by Eq.(25). This is because the cylindrical CMC slices act as barriers to the non-cylindrical CMC slices. 

As C approaches C+, each of the two CMC slices will develop long cylinrical regions. The one from null infinity will 
run along, but just above the surface with r = Rf while the one from r — will run just below. The closer to the 
critical value, the longer the cylinders. 

When C = C+ we get a sudden change. Instead of having two solutions with long cylinders we have five. Two come 
from the left and right null infinity, respectively, and asymptote (from above) to infinite cylinders of radius r = Rf . 
Two others come, left and right, from r = and asymptote from below to the same cylinders. The fifth solution is 
the r = Rf cylinder itself. 

When C exceeds the critical value we get another change. The polynomial becomes everywhere positive. This 
means that the CMC slice cannot have any extremum. It must run all the way from r = to r = oo. If C > we 
will have two CMC slices, one from the left future null infinity which runs into the future singularity and a mirror 
one from the right future null infinity. 

Starting from the umbilical slice, holding K fixed, and let C become negative, we get the opposite behaviour. 
The slice from null infinity to null infinity moves backwards in time, while a new CMC slice emerges from the past 
singularity and goes back to it. As C approaches a negative critical value C~ , the two roots of the polynomial 
approach one another and coincide at a radius R~ < 3m/2. This is the smaller root of Eq.(22) and C~ is the value 
of C given by Eq.(26). 

We conjecture that the slicings we have described for fixed K and for C in the range C~ < C < C+ form three 
separate foliations, one for r < R~ near the past singularity, one for r < i?+ near the future singularity, and third 
foliation formed by the slices that run from one null infinity to the other. We further conjecture that these three 
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foliations cover the entire extended spacetime. 



V. DIFFERENTIATING THE HEIGHT FUNCTION 



Let us consider the foliation that runs from null infinity to null infinity. Each slice has the same value of K but C 
spans an interval. We could use the value of C as a label on the slices but we want to use some time coordinate as 
label. The obvious choice is the 'time at infinity'. This is given by 

roc IKr _ C_\ j roc - c) dr 

t(C) = / 7 r2)dr = / } 3 > (41) 



(1-2*) y/l-2? + (Ig-$) 2 ^ (l- 2 ?) ^r±-2mr* + (Zf-Cy 

This has three divergences. The first is due to the (1 — 2m/r), which diverges at the horizon. This can be integrated 
through in the Cauchy principal value sense. A similar divergence arose in [6] . The second is due to the fact that the 
polynomial inside the square root vanishes at r — r m i. This is the definition of r m i- This is not a problem because 
the polynomial goes to zero linearly at r m ;. Therefore the integral is of the form J dx/^fx, which is regular at x = 0. 

The third divergence is due to the fact that the integral itself diverges as r — > oo. This has to be because the slice 
goes to null infinity. To leading order the integral becomes 



fdr / „ 9 . 



which is just the flat spacetime mass hyperboloid. If we want a finite time label on the CMC slices the obvious 
thing to do would be to subtract off the leading flat space divergent expression. Unfortunately, the difference still 
logarithmically diverges (like 2m In r) . If we want a finite expression it is better to subtract off the height function of 
some favored slice of the foliation itself. One obvious choice is to pick the umbilical slice (the C = 0) slice. Therefore 
a natural time label is given by 

roo - C) dr r oo Ellrlr 

r(C)= } 3 I = / 3 (43) 

Jr ml (!_ M) ^ r 4-2mr3 + (Zf-C) 2 Jr - (l - 3s) yV 4 - 2mr^ + (^ff 

This, from the argument given above, is finite for all C < C*. As C — > C*, we have that r m ; — > R*. At this point, 
the two roots of the polynomial coincide, and the slope of the tangent to the polynomial at r = r m ; goes to zero. The 
integral close to r m i approximates J adx/y/sx where a is some constant and s is the slope. Integrating this over some 
small but finite fixed interval (0, Ax) we get a contribution to r of as/ Ax/2^/s. As s — > 0, this contribution becomes 
unboundedly large. Therefore we get 'collapse of the lapse' in the interior. The foliation moves only a finite distance 
at the center to reach r = i?* while the passage of 'time at infinity' becomes unboundedly large. 

At the critical point both \/h' and the first derivative ofl/h' vanish at the throat. The coefficient in the exponential 
decay is nothing more than the second derivative of 1/h! at the critical point. This is the dominant term in any 
expansion of the time function near the critical point. The rest of this article is devoted to demonstrating this. 

We wish to investigate the behaviour of the central lapse. In [6] we discussed the situation where we had a foliation 
defined by some time function r with lapse a. Say we are given a vector field The projection of £ normal to the 
time slice (call it N) is given by 

N = a^V => a = N {CV ^y 1 . (44) 
If we choose £ to be the Killing vector, we know what is from Eq.(29) and we also can write 



n-1 ( dr\ 1 dh ( dr\ 1 

rv - r) ={dc) dc r ^ a= {dc) 



dr\ , T dh 
N dC 



(45) 



To evaluate expression (45) we need to differentiate the height function with respect to C. This looks to be highly 
unpleasant. The square root in the denominator is promoted to 3/2 power so the integral has a term dx/x 3 ^ 2 which 
diverges at the origin. Further, r TO ; depends on C so there will also be an end-point variation. This will take the 
integrand (which is infinite) outside the integral sign. We know dr/dC must be finite so these infinities must cancel. 
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A very similar problem arose in [6] and a way was found around it. This essentially involved an integration by parts 
before differentiating and much more malleable expressions were found. We can repeat this trick. 
We begin by defining the following function 



r 4 -2mr 3 + 



dr 



This is constructed so that dJ/dC = h. Now rewrite J as 



(46) 



2 r 
~3j WZ 



r 4 - 2mr 3 + - c)' 



dr 



(l - 3m) (4r 3 - 6mr 2 + 2Kr 2 [^f - C] ) 
This now can be integrated by parts to give 



(47) 



J 



3 (1 - 3s) (4r 3 - 6mr 2 + 2Kr 2 [Zf - C] ) 
Kr 3 



2mr J 



C* 



(l _ 3s) (4r 3 - 6mr 2 + 2i^r 2 [*j± - C] ) 



dr. 



(48) 



We need to differentiate this twice with respect to C to get dh/dC. We will do this in two parts. Let us call the 
not-integral part J\ and the integral J 2- 
We get 



dJi 
~dC 



= 2 



r 4 - 2mr 3 + - c)1 * - C 



2mr 3 + - C) 



(1 - M) ^ 4r 3 _ g mr 2 + 2Xr 2 - C] ) 

2" 



ifr 2 



3 (! _ 3m) ( 4r 3 _ 6mr 2 + 2 £>2 _ C ] )' 



(49) 



d 2 Ji 
dC 2 



= - 2 



- 2 



r 4 -2mr 3 + (i^-C) 2 * - C 



'I _ 3s) (4r 3 - 6mr 2 + 2ifr 2 - C]) 

1 



+ ■ 



(l - 2a) (4r 3 - 6mr 2 + 2Xr 2 - C] ) 

1 

r 4 - 2mr 3 + - df\ ' Kr 2 - C) 



(1 _ 3s) ( 4r 3 _ 6mr 2 + 2 #>2 _ ) 





r 4 -2mr 3 + (^-C) 2 


2 ^v 


16 







(l - 2m) ^ 4r 3 _ g mr 2 + 2X r 2 _ C] ) 



,3' 



(50) 



One interesting property of d 2 J\/dC 2 is that it vanishes for large r. This means that it does not contribute to dr/dC. 
Note also that the first term in d 2 J\/dC 2 diverges as r — > r m ;. However, we must remember that to compute a we 
multiply by N which goes to zero in the matching fashion so that everything is regular. Further, only the first term 
is finite at the throat, all the other ones vanish. 
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Now we can work out 



dC 



(1 - 221) ( 4r 3 _ 6mr 2 + 2 Kr 2 [Kf-C]) 



(51) 



From the definition of r TO / as the zero of the polynomial, Eq.(38), it is clear that 

(^f-°) 2 =( 2mr3 - r4 )r ml ->(nr- C ) =-y/(2rnr>-r*) rmi . (52) 
^From Eq.(21) we have that K = (2i?» — 3m) / y / 2mRl — i? 4 . When these are substituted into Eq.(51) we get 

1 



dC 



2r-3m-(2 J R,-3m)y^PJ 
A natural variable to use (as in [6] ) is 5 = r m l — R* ■ We then get 



(53) 



N 



dh_ 

dC 



1 



2mi?„ - R 2 



(3m-2fi.) 2 
2mR,-Ifj 



s (2R 2 -8mR*+9m 2 )5' 



(54) 



The polynomial 



D = 4r 3 - 6mr 2 + 2Kr 2 



Kr 3 



C 



(55) 



in Eq.(51) is the first derivative of the sextic polynomial of Eq.(38). In general it does not vanish at r TO ;. However, 
we can see that 



D&2R 2 



2 + 



(3m- 2R„) 2 



2mR* - Rl 



(56) 



and thus, as expected, goes to zero as C — > C* . 

Now we need to look at the integral part of J as this is what gives us dr/dC. 



J 2 = 



r 4 - 2mr 3 + ( - C 



-3r 2 + 7mr - - 3m 2 + imI £ r * + (Kr - mK) C 

[2r 3 - 7mr 2 + ^f- + 6m 2 r - + (2mKr - Kr 2 ) C] 



7 dr (57) 



dj 2 
dC 



= -2/ 



r 4 - 



\ 3 / J V 3 / [2r 3 -7mr 2 + -gfi+6m 2 r- 2m f ^ + (2mKr-Kr 2 )c] 

3 



■ dr 



r 4 — 2mr 3 



(Kr-mK) 



[2r 3 -7mr 2 + S^±+ern 2 r- 2 ™tf^ + (2mKr-Kr 2 )C'] 



rdr 



-II 



r 4 — 2mr 3 



( Ell 
V 3 



C 



(-3r 2 +7mr- hK \ Tl -Zm 2 + imI f^ +(Kr-mK)C*) (2mKr-Kr 2 ) 
[2r 3 -7mr 2 +-^+6m 2 r- 2m f r4 +(2miCr-Kr 2 )c] 3 ' 



(58) 



d 2 J 2 _ dr _ ' 
dC 2 ~ dC~ + 



r 4 - 2mr 3 + I ^ C 



-C) x 



-3r 2 + 7mr - - 3m 2 + 4mJ f r " + (iiTr - mJf) C 



+ eight other terms 



[2r 3 — Imr 



2^Ell + Q m 2 r - 



(2mKr - Kr 2 )C\ 



■dr. 



(59) 
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All the nine terms in Ew.(59) fall off like 1/r 3 . Therefore each of these terms is finite. We also know that 
dr/dC — > oo as C — > C*. The term we have isolated is the term which generates this behaviour. All the other terms 
remain finite. To estimate the blowup we need to understand the behaviour of it near r TO ;. It is useful to shift the 
origin of coordinates to r = r cr u = Therefore we define y = r — R*. We know K = (2i?# — 3m)/ \/2rnR% — 
and we also write C = C* — e, where C* = (3mRl — R^)/3y / 2mRl — R%. 

We now write out the polynomial Eq.(38) in terms of (R*, m, y, z) to give 

4 2R 3 -8mR 2 + 9m 2 R* 2 12i? 2 _ 56i?„m + 54m 2 , 17.R 2 - 54™^ + 45m 2 , 

r = — v H v H 1 v 

2m- R* y 3(2m-R*) 3a(2m-i?*) y 

AR 2 - \2R,m + 9m 2 r - fil 2i? 3 - AR 2 m , 4i?» - 6m r , , , n 

+ •»«.(»»-«.) ' 6R -" 5 + " 8 i + jsseh? + ' + ^r.-u^ -" + «•*" + " le - < 60) 

The polynomial begins at y 2 because we know that when e = both the polynomial itself, and its first derivative 
vanishes at y = 0. More generally, we know that the polynomial vanishes when r = r m /, i.e., when y — r m i 
If e is small, if we are close to the critical value, then 

2Rl - 8mi? 2 + 9m 2 i?* r2 2Rl - AR 2 m 2Rj - 8mi?J + 9m 2 i? 2 r2 , x 

— *— S 2 + * e ~ => e pa i =(5 2 . (61) 

2m -R* v /2mi?» - i? 2 (4# 2 m - 2Rl)^j2mR* - R 2 

Further, near r = r m ;, we find that the polynomial approximates 

4 2i?g - 8mR 2 + 9m 2 i?, 2 2 

r - = — 2^: — (62) 

The poynomial Z) of Eq. (55) is the first derivative of Eq. (60) . Thus we have 

dD = 2Rl-8,„Rl + <>„,:> R 

dr 2m -R* yy ' v ' 

We also need to approximate the other terms in Eq.(59). The denominator equals (1 — 2m/r) 2 D 2 /A. When we write 
this out in terms of (i?* , m, y) we get 

2r3 . . . « ( 1 - ^V 2 ^- 8 ^ + 9 ^V + 0(y 3) (65 ) 



and 



The polynomial in the numerator of Eq.(59) is 



ld_ 
2~dr 



l - 2 -y) D 



(66) 



Therefore 



2 , 1„ 2m N 2i?2-8mi? 2 + 9m 2 i?» 



-3r 2 + = --(1 - -) * 2m : R + 0(y). (67) 



Therefore, the term in Eq.(59) which blows up as 5 — > is dominated by 

d 2 J l= d Lr r {^- C ) ( 2m -R* \ * 1 

dC 2 dC~ + J S 2 -f-\ \2R\ - 8mRl + 9m 2 R J y^^y^~P V ' { ' 

We know that 

V 2 SL _ 1 J * Ri (69) 
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and 



J y 2 V^^ V v s2 ' 

Therefore we have 

d 2 J 2 _ dr_ _ / 2m -R* \ 1 i£ 



VI. COLLAPSE OF THE LAPSE 

We now have calculated the various terms necessary to compute the central lapse. ^From Eq.(54) we have 



dC 



2mR* - Rl 



[2R 2 - 8mR sr +9m 2 )5' 



^From Eq.(61) we have 



e = C. - C ■. 



2Rl - 8mR* + 9m 2 



(4m- 2R*) y j2mR st - R 2 



-J 2 . 



We also have, from Eq.(71) 



dT 



2m- R* 



dC \ 2RI - 8mR 2 + 9m 2 i?» J 5 2 ' 



Rt 



We can differentiate Eq.(73) to get 



dC 



2Rl - 8mR* + 9m 2 



dS (2m- R*)y/2ma - R 2 



S. 



We multiply Eq.(74) by Eq.(75) to get 



dr dC dr 
dC~dJ = d8 



Rt 



2R 2 - 8mi?* + 9m 2 / 5 



Integrating Eq.(76) gives 



r = 



Rt 



2R 2 - 8mi?* + 9m 2 



In 5 + A 



where A is a constant, or 



S = exp 



2R 2 - 8mi?„ + 9m 
RT~ 



2\ 2 



(r-A) 



^From Eq.(45) we have 



dr\ 1 ^dh 

: \dc) N dc 



Therefore to compute the central lapse we need to divide Eq.(72) by Eq.(74) to get 



/ dr \ \ T dh 

a =[dc) N dc 



2m — i?» 



2RI - 8mR 2 + 9m 2 i?* V 2m - R 



2Rt - 8mRi + 9m z R 



2D \ 2 



R 2 



2Rl - 8mR„ + 9m 
2mi?3 _ R i 



2\ 2 



12 



On substituting in Eq.(78) we get 



a = B cxp 



2i? 2 - 8mR* + 9m : 



where B is a constant which equals 



B = 



2Rl - 8mR* + 9m 2 
2mRl - Ri 



exp 



2Rl - 8mR* + 9m 2 
Ri 



A 



(81) 



(82) 



We have not evaluated A so we cannot compute B. 

It is clear from Eq.(78) that A sets the zero of r. When we compute the collapse of the lapse for the maximal case 
the moment of time symmetry slice sets a natural zero for the time function. In the CMC case we discuss here we 
cannot use the Killing time at infinity because it is infinite. As discussed in the beginning of Section V we have to 
normalize it by setting the zero of time to be that of the umbilical slice, the slice which has C = 0. In some sense, 
this is the analogue of the moment of time symmetry slice, but at the same time it is somewhat arbitrary. This 
arbitrariness will be reflected in the constant A. It is also quite difficult to compute. We would need to evaluate 
integral (43) as we approach the critical slice. The leading term should agree with Eq.(77) but we also need to compute 
the next term, which will give us A. 

This indicates a different way of computing the exponent. Let us look at h' as given by Eq.(35), or rather, let us 
look at 1/h' 2 near the critical point: 



1 

V 2 ~ 

When this is substituted into Eq.(43) we get 

r(S) 

We have 



2Rl - 8mR* + 9m 2 , _, 



Ri 



Rl 



dy 



s ^2Rl - 8mR* + 9m 2 sjy 2 - S 2 



I 



dy 



In 



y + \Jy 2 -6 2 



In 6. 



(83) 



(84) 



(85) 



Therefore we reproduce Eq.(77). 

There are many ways of deriving the critical exponent. In a companion paper, [11], we offer a very different 
derivation, based on an explicit formula for the lapse function for spherical CMC slices. The computation given here, 
and the calculation there entirely agree. 

The calculation here is closely modelled on the computation of the critical exponent for the maximal foliation given 
in [6]. This allows us to perform a number of internal consistency checks at various points in the calculation by 
reducing to the K — situation. We get agreement at each stage. This agreement is not trivial because the key 
integration by parts to obtain an explicitly finite derivative of J differs in the two cases. 



VII. FOLIATIONS WITH NON-CONSTANT K? 



As is clear from the discussion, one has three parameters to play with in constructing spherical CMC foliations or 
slicings in the Schwarzschild solution. One can change K, one can change C, one can drag a slice along the 'timclike' 
Killing vector. Of course, ons can change more than one of these at once. The foliation we have focused on is one 
where we kept K fixed, changed C and eliminated the Killing freedom by considering the slices where the minimal 
surface coincided with the t = axis in standard coordinates. 

We could consider any one of these slices and drag it along the Killing vector. For the slices which run from null 
infinity to null infinity, one with \C\ small, one would get a slicing which looks somewhat like the standard t — 
constant slicing of Schwarzschild. It would rise up along one null infinity and sink down on the other while the throat 
ran along one of the R = constant lines in either the upper or lower quadrant. These would not form a foliation, each 
slice crosses each other in the interior, the lapse function will vanish on the throat and be positive on one side and 
negative on the other. If we drag one of the slices which run from null infinity to the singularity, one with |C| large, 
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along the Killing vector we get a foliation, the slices do not cross, but each one ends on the R = singularity. Such a 
foliation would cover the upper half of the right hand quadrant and all the upper quadrant if one picked K > and 
C > 0. Other patches could be covered by choosing other options for K and C. 

In closed cosmologies, on the other hand, we are used to CMC foliations where the value of K changes. If the 
cosmology goes from a big bang to a big crunch we might expect to have a foliation which goes from K = +oo, at 
the big bang, through K — 0, the moment of maximum expansion, and monotonically to K = — oo, the big crunch. 
It can be shown that no such foliations with varying K exists in the Schwarzschild solution. 

A first try would be to consider the slicing where one changes K while keeping C fixed. Such slices always cross 
each other. Consider two slices, one with K = 0, C = (this is the standard moment-of-time-symmetry, t — 0, slice 
through Schwarzschild), and the other with K = 1, C = 0. As we discussed previously, this slice starts at future null 
infinity, crosses in the middle below the bifurcation point and rises up again to the other future null infinity. This 
slice crosses the first slice twice. Choosing a different, fixed value of C will not change this behaviour. Therefore, if 
we want a foliation with varying K, we need a nonconstant C. 

One almost has such a foliation. Consider the slicing where one changes K while simultaneously changing C such 
that C = 8m 3 K/3. Each of these slices has its throat at R — 2m. Each of these slices runs through the bifurcation 
point and so they must all touch there. They do not cross, however, and this is their only point of contact. This 
slicing covers all of the left and right quadrants. One might hope that by slightly changing C, one could spread the 
slices apart along the vertical t — axis and thus convert this slicing into a foliation. This cannot be done, as we 
show below, if we want to allow K to be unboundedly large. 

The standard work on the way CMC slices act as barriers was written by Brill and Flaherty, [12]. Among other 
results, they show that two slices with the same value of K cannot touch at a single point. Further, if two CMC 
slices do touch at one point, the slice to the future must have a larger value of K than the other one. This restriction 
strongly restricts the behaviour of CMC slices in the Schwarzschild solution. 

Let us assume that we have a CMC foliation of a Schwarzschild solution which starts off at the moment of time 
symmetry slice and moves up. Consider one slice of this foliation. This slice will have positive K — Ks and 
C > 8m 3 K/3. This slice will have a throat with some radius Rs- The cylindrical slice with this given radius is also 
a CMC slice (with K = K\, say) and touches the other CMC slice at one point, the throat. Therefore we must have 
that Ks < K\. Since we assume we have a foliation, we must have that Rs monotonically decreases and hence also 
K\. It passes through zero at Rs = 3m/2. However, we expect that Ks to be increasing as the foliation moves 
forward in time so we will eventually run into a situation where K s — Ki, which Brill and Flaherty forbid. Therefore 
any spherical CMC foliation of Schwarzschild cannot have unboundedly large values of K. 
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